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TITLE:  tAertxcal  Structure  of  Flow  Through  Tarifa  Narrows 


The  research  undertaken  in  this  project  is  part  of  a  larger  effort,  in 
collaboration  with  N.  Pettigrew  (Uhiversity  of  Maine,  Orono) ,  in  which  the  PI 
seeks  a  description  and  a  dynamical  explanation  of  the  steady,  tidal,  and 
internal  wave  components  of  the  current  in  the  Tarifa  Narrows  of  the  Strait  of 
Gibraltar.  With  this  effort,  the  PI  continues  his  summer  research  program 
while  teaching  undergraduates  fulltime  during  the  academic  year  at  Norwich 
Uhiversity.  The  research  funded  by  this  CXR  grant  permitted  the  PI  to  work 
the  summer  of  1992  at  the  Ocean  Process  Analysis  Laboratory  (Uhiversity  of  New 
Hampshire;  W.  S.  Brown,  director).  The  PI  was  involved  primarily  with  the 
analysis  of  Acoustic  Dopier  Current  Profiler  (ADC]P)  data  taken  by  Pettigrew 
and  collaborators  (Proc.  of  Ocean  ’87^  1987)  during  the  Gibraltar  Experiment 
(GIEEX)  (Kinder  &  Bryden,  BOB  Transact,  ABU,  1987). 

The  format  for  this  final  technical  report  will  outline  the  principal 
objectives  of  the  research  and  describe  the  progress  made  on  each.  In 
general,  the  research  is  not  finished  and  continues  to  occupy  the  attention  of 
1  the  PI.  In  particular,  the  PI  will  work  to  finish  several  of  the  uncompleted 
objectives  during  future  summers. 

I 

The  PI  is  most  appreciative  of  the  funding  provided  as  he  was  able  to 
advance  on  many  aspects  of  the  project. 


GEUEZ777VE'  Is  a  refinement  of  the  response  method  of  tidal  analysis  for  use 
with  current  records  strongly  contaminated  by  internal  bores. 
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The  unique  ADCP  data  set  from  Tarifa  is  profoundly  characterized  by  the 
passage  of  a  tremendous  internal  bore  generated  semidiurnal ly  at  the  nearby 
Camarinal  Sill.  To  avoid  contamination  by  the  bore,  a  modification  of  the 
response  method  of  tidal  analysis  has  been  developed  to  extract  the  steady 
component  of  flow  and  the  barotropic  tide  from  the  dynamically  complex 
current;  the  bore  is  left  as  the  residual  of  the  analysis.  This  decomposition 
permits  the  steady  component  and  the  internal  bore  to  be  examined  as  isolated 
dynamical  entities. 

The  development  of  this  technique  \nas  finalized  by  the  PI  during  the 
summer  of  1992  while  funded  by  this  grant.  A  photocopy  of  an  article  which 
finally  appeared  (Hyde  &  Pettigrew,  JGR  99,  #Cir,  1993)  is  attached. 

Figure  1  (taken  from  Hyde  &  Pettigrew,  1993,  page  22,755)  shows  an 
example  of  difference  between  the  standard  response  method  and  the  modified 
response  method  of  tidal  analysis.  The  two  methods  produce  drastically 
different  values  for  the  amplitude  of  the  tide,  the  phase  of  the  tide,  and  the 
spend  of  the  mean  flow  as  well  as  the  structure  of  the  internal  bore,  and 
therefore,  quite  different  dynamical  views  of  the  Strait. 


OBJECTH/E  2:  a  completion  of  the  decomposition  of  Tarifa  ADCP  time  series 
into  their  steady,  tidal,  and  internal  bore  components. 

Approximately  257.  of  the  Tarifa  ADCP  data  was  analyzed  confidently 
during  the  summer  of  1992  under  the  funding  of  this  grant.  In  particular,  the 
2  periods  of  spring  tides  at  the  deep  Tarifa  instrument  were  finished. 

However,  the  neap  tides  are  much  more  F)roblematic  and  require  greater 
dynamical  insight  before  they  can  be  decomposed  with  confidence;  the  same  is 
true  of  the  entire  ADCP  record  in  the  shallows  by  Tarifa  Harbor. 

Nevertheless,  this  analysis  has  provided  a  very  interesting  data  set  showing 
the  vertical  variation  of  the  mean  flow  and  the  internal  bore  during  the 
periods  of  spring  tides. 

Figures  2  and  3  show  the  results  of  the  analysis  during  a  period  of  ten 
semidiurnal  tidal  cycles  of  the  spring  tide  in  early  April,  1986  in  the  deep 
Tarifa  Narrows.  Figure  2  shows  a  time  series  of  the  variation  with  depth  of 
the  velocity  associated  with  a  composite  of  the  ten  internal  bores  which 
passed  the  instrument  during  the  period.  Figure  3  displays  the  vertical 
variation  of  the  north-south  and  east-west  components  of  the  steady  flow 
during  the  period. 


CR7£CT/VE  J:  an  estimate  of  the  transport  of  Atlantic  water  into  the 
Mediterranean ; 

An  estimate  by  Pettigrew  &  Hyde  (Phys  Otean  Ssa  Strait^  NATO/ASI,  1990) 
of  the  transpiort  of  Atlantic  water  into  the  Mediterranean  showed  that  the 
contribution  to  the  transport  by  the  internal  bore  is  comparable  to  that  of 
the  steady  flow. 

The  funding  of  this  grant  permitted  a  considerable  improvement  in  the 
earlier  estimates.  Applying  the  modified  response  method  to  the  upper  layer 
of  inflowing  Atlantic  water  to  isolate  the  bore,  the  steady  flow,  and  the 
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FIGLRE  1  Omparison  of  the  standard  response  method 
of  tidal  analysis  (A  &  B)  with  the  modified  response 
method  of  tidal  analysis  (D  &  E).  Tarifa  DAPCM  data 
at  40  meters  depth. 
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tide,  the  mass  transport  of  Atlantic  water  through  the  Narrows  was  carefully 
decomposed.  When  averaged  over  a  spring-neap  cycle,  it  was  foind  that  the 
transport  associated  with  the  steady  flow  is  0.56  Sverdrups  while  that 
associated  with  the  bore  is  0.43  Sverdrups.  The  barotropic  tide's 
contribution  is  negligible.  Thus,  the  importance  of  the  bore  noted  in  the 
earlier  investigation  was  reaffirmed.  It  was  also  found  that  the  transport 
during  a  neap  tide  (1.06  Sverdrups)  is  slightly  greater  than  during  a  spring 
tide  (0.93  Sverdrups). 

The  PI  has  completed  the  computations  required  for  this  objective.  A 
manuscript  (with  N.  Pettigrew  as  lead  author)  now  needs  to  be  prepared. 


OBJECTI\^  4i  a  dynamical  model  of  the  observed  vertical  structure  of  the 
steady  component  of  the  flow  at  Tarifa; 

Figure  3  shows  the  depth  dependency  of  the  eastward  (triangles)  and  the 
northward  (crosses)  components  of  the  steady  flow  during  a  period  of  10  spring 
semidiurnal  tidal  cycles.  There  is  considerable  shear  and  turning  in  the  mean 
flow  across  the  Atlantic-Mediterranean  interface  suggesting  the  existence  of 
boundary  layers  straddling  the  interface. 

At  the  suggestion  of  Pettigrew,  the  PI  investigated  the  vertical 
structure  of  the  mean  flow  using  Ekman  boundary  layer  theory.  The  PI  fitted 
the  data,  as  shown  in  Figure  3,  with  a  simple  Ekman  model  which  balances 
pressure  gradient,  Coriolis,  and  eddy  viscous  forces  in  the  interior  of  the 
fluid  layers  and  imposes  appropriate  boundary  conditions  on  the  velocities  and 
fluid  stresses  at  the  inter-face,  at  the  sea  surface,  and  at  great  depth.  The 
fitting  parameters  are  the  wind  stress,  the  pressure  gradients  in  the  two 
layers,  the  depth  of  the  interface,  and  the  mixing  coefficient.  The  fit 
(Figure  3)  is  good  and  the  values  of  the  fitting  parametei-s  are  in  reasonable 
agreement  with  the  physical  reality  of  Tarifa.  This  success  provides 
encouragement  to  continue  this  analysis. 

Roughly  5051  of  the  work  required  for  this  objective  was  accomplished  by 
the  PI  during  the  summer  of  1992  with  funding  from  the  CNR  grant.  There  are 
details  to  be  considered  and  there  is  an  analysis  of  the  data  from  Tarifa 
Harbor  to  be  attempted  before  a  manuscript  can  be  prepared  (with  Pettigrew  as 
lead  author). 


£BJEC7J(C  5:  a  dynamical  model  of  the  temporal  evolution  and  vertical 
structure  of  the  internal  bore  at  Tarifa. 

The  composite  of  the  horizontal  velocity  field  associated  with  the 
internal  bore,  shown  in  Figure  2,  is  broadly  consistent  with  the  known 
characteristics  of  a  tuor-layer  solitary  wave,  namely,  longitudinal  motion  in 
the  direction  of  propagation  (eastward)  in  the  thinner  upper  layer  and 
opposite  the  direction  of  propagation  in  the  thicker  lower  layer.  However, 
there  are  also  considerable,  complex  transverse  motions.  In  particular, 
during  the  latter  part  of  the  passage  of  the  wave,  it  appears  that  the 
tranvtmcrse  velocity  is  primarily  to  left  Cright]  of  the  propagation 
direction  in  the  lower  [upper]  layer. 
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TIME  (HOURS) 


FlGLffE  2  Time  series  composite  (average  of  10 
semidiurnal  tidal  cycles)  of  the  horizontal  velocity 
associated  with  the  internal  bore  as  it  passes  Tarifa 
Eastward  current  to  right;  northward  current  up.  Tim 
measured  from  the  zero-crossing  of  the  tidal  pressure 
record  at  Tarifa.  Ghset  of  the  bore  at  about  4  hours, 
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FIGUFE  3  NorthMard  and  eastward  components  of  the 
steady  flow  during  a  Tarifa  spring  tide. 


An  attempt  to  explain  the  dynamics  of  the  transverse  motions  of  the  bore 
has  not  boon  successful  so  far.  An  asymptotic  analysis  of  solitary  waves 
(weakly  dispersive,  weakly  nonlinear,  long  waves)  influenced  by  rotation, 
which  the  PI  submitted  to  JFM,  was  deemed  inafspropriate  by  reviewers. 

Instead,  the  modified  Kadomtsev-Petviashvilli  (MKP)  equation  (Grimshaw  & 
Melville,  Stud  Appl  Math^  1989)  was  suggested.  However,  it  is  not  immediately 
obvious  that  the  behavior  in  Figure  2  is  explained  by  MKP  theory. 

Considerably  more  theoretical  work  needs  to  be  done  in  this  regard.  Perhaps 
due  to  its  evolution,  the  bore  exhibits  several  velocity  fronts  and  a 
complicated  depression  of  the  interface  which  will  frustrate  any  definitive 
explanation  of  the  bore's  dynamics. 

At  a  minimum,  the  PI  will  attempt  to  publish  the  data  showing  the 
structure  of  the  transverse  motions  of  the  bore. 
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A  Modification  of  the  Response  Method  of  Tidal  Analysis 

Richard  A.  Hyde 

Pkfsiu  Dtpartmua.NorwUk  Umvtrtiiy,  Norlhfield,  Vermont 
NEALR.PETnCREW 

Departmtnt  of  Oetonotropky,  University  ef  Maine  at  Orono 

An  easily  anplemeined  exieatiaa  of  the  standard  response  method  of  tidal  analysis  is  outlined.  The 
modificatian  improves  the  ezuactioii  of  both  the  steady  and  the  tidal  compooenu  from  proUemaiic  time 
setiet  by  calculating  tidal  respaose  weights  uncontaminated  by  misting  or  anomalous  data.  Examples  of 
time  series  eontaining  data  gaps  and  anomalous  events  are  analyzed  to  demonstrate  the  applicability  and 
advantage  of  the  proposed  metlMd. 


INIKCXXlCnON 

The  analysis  of  tidal  cunenu  may  be  confounded  by  an 
inability  to  separate  clearly  the  surface  (barotropic)  tide  and 
the  internal  (baroclinic)  tide.  This  is  particularly  true  in 
regions  of  internal  tidal  generation  such  as  the  Strait  of 
Gibraltar.  Data  obtained  during  the  Gibraltar  Experiment 
[Kinder  and  Brjden,  1987]  include  enormous  velocity 
surges  associated  with  the  passage  of  a  spectacular  internal 
bore  [Pettigrew  and  Hyde,  1990]  that  is  approximately 
phase  locked  to  the  surface  tide.  Surge  velociUes  often  are  of 
the  order  of  a  meter  per  second  and  may  substantially  exceed  all 
of  the  other  flow  componenu  near  the  pycnocline.  In  a  fash¬ 
ion  'common  to  the  generation  of  most  baroclinic  tides 
[Rattray,  I960:  Wunsch,  1975],  this  large  bore  is  produced 
u  a  consequence  of  the  surface  tide  driving  a  sharply  stratifted 
water  column  over  the  Camaiinal  silL 

Intermittency  of  anqrlitude  and  phase  [tPenscA.  1975],  due 
presumably  to  temporal  and  spatial  variations  in  hydiogrqdiic 
conditioiu  either  in  the  generation  area  or  in  the  radiation 
field,  often  characterizes  internal  tide  oscillations.  For  very 
long  current  records  far  from  the  generation  region,  the  baro- 
clii^  component  of  the  tide  often  appears  as  little  irwre  than 
iwise  at  the  barotropic  tidal  frequeiKies  and  may  therefore  be 
removed  with  a  standard  response  analysis  [Munk  and  Cart¬ 
wright,  1966]  of  the  tide.  At  the  other  extreme,  in  regions 
where  local  internal  tidal  generation  is  large,  the  cutrenu  in 
the  tidal  band  may  have  a  very  complex  structure  that  is  phase 
locked  with  the  surface  tide.  The  case  of  the  Strait  of  Gibraltar 
is  a  particularly  dramatic  example  of  this  latter  extreme. 
There,  a  very  large  artqrlitude  bore  severely  distorts  the  veloc¬ 
ity  signature.  (See  the  time  series  rimwn  in  Hgure  Id  which 
be  discussed  thoroughly  in  section  3.)  Standard  tidal 
analyses  therefore  produce  unreliable  amplitudes  and  jdiases. 
Moreover,  the  baroclinic  residual  of  the  analysis  bears  little 
resemUance  to  known  characteristics  of  the  bore  itself.  Thus 
indiscriminate  use  of  the  traditional  response  mediod  (or  of  the 
harmonic  method)  of  tidal  analysis  may  lead  to  fallacious 
views  of  the  dynamics  of  the  strait  A  modifted  qtproach  is 
needed. 

The  response  method  of  tidal  analysis  was  introduced  [Munk 
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and  Cartwright,  1966]  as  an  alternative  to  the  harmonic 
method  of  tidal  analysis  [Darwin,  1883].  This  method,  devel¬ 
oped  further  in  the  work  of  Cartwright  et  al.  [1969]  and 
Zetler  and  Munk  [1975],  is  based  upon  the  relationship 
between  an  input  7(1)  and  the  output  v(/)  of  a  linear  system. 
Namely, 

y(t)  =  W(0nt-^dT  (1) 

where  1F(t)  is  an  impulse  response  function.  The  Fourier 
transform  of  (1)  results  in 

K®) -H(o>)7Xo>X  (2) 

whae  U  is  variously  referred  to  as  the  transfer  function  or  the 
admittance  and  is  obviously  the  Fourier  transform  of  the 
coherent  ratio  of  output  v  to  input  T. 

The  response  method  offers  several  advantages  over  the 
harmonic  method  [Zeller  et  al„  1979;  Zetler,  1987],  includ¬ 
ing  greater  accuracy,  clearer  separation  of  solar  gravitational 
and  solar  radiational  effects,  and  better  statistical  stability. 
The  last  feature  is  perhaps  the  greatest  distinction  between  the 
two  methods.  While  a  reliable  estimate  of  tidal  constituents 
by  harmonic  (or  spectral)  analysis  requires  long-term  observa¬ 
tions.  the  response  method  can  analyze  a  relatively  short 
record.  It  achieves  both  stability  and  resolution  by  utilizing  a 
long-term  reference  series  upon  which  high-resolution  analy¬ 
sis  can  be  performed  and  then  calculating  a  smooth  traiufer 
function  [Munk  and  Cartwright,  1966]  linking  the  reference 
series  to  the  observed  time  series.  (See  May  [1979]  for  a  brief 
but  lucid  discussion.)  As  a  resulL  the  response  method,  al¬ 
though  not  universally  employed,  has  become  the  standard  for 
the  analysis  and  prediction  of  tidal  currents. 

This  document  presenu  a  modification  of  the  response 
method  of  tidal  analysis  which  seems  quite  useful  for  isolating 
anomalous  events  (such  u  passages  of  a  bore)  and  preventing 
them  from  contaminating  tidal  analyses.  The  procc^e  is  an 
extension  of  and  improvement  upon  the  preliminary  work  of 
Pettigrew  and  Hyde  [1990]  for  discriminating  between  the 
Gibraltar  internal  bore  and  iu  engendering  tidal  current. 
Termed  "the  selective  response  method  of  tidal  analysis,”  it 
alters  the  response  method  formalism  so  that  only  those  data 
that  are  clearly  uncontaminated  by  anomalous  eventt  and  are 
clearly  representative  of  the  steady  flow  and  the  tidal  oscilla¬ 
tion  are  included  in  the  emnputation  of  the  traiufer  function. 
While  development  of  the  selective  respoiue  method  was 
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motivated  by  the  diHicuItiet  associated  with  the  Gibraltar 
internal  tide,  the  method  lends  itself  naturally  to  the  analysis 
of  velocity  records  with  missing  data.  Because  such  records  are 
so  common  and  because  preliroinaiy  tests  have  been  so  suc¬ 
cessful.  it  may  well  be  in  its  application  to  data  with  gaps  that 
the  modified  method  will  find  its  greatest  utility  to  ocean¬ 
ographers. 

SELBCnVE  RESPONSE  MEIHQD  OF  UDAL  ANALYSIS 

The  selective-response  method  of  tidal  analysis  is  founded 
upon  the  following  two  modifications  of  the  stmidard  response 
method:  (1)  an  alteration  of  the  calculation  of  the  response 
weights  in  which  only  selected  data  determine  the  optimal 
(least  squares)  fit  of  the  tide  and  the  mean  and  (2)  an  adjustment 
in  the  determination  of  the  appropriate  number  of  respoiue 
weights. 

The  selective  response  method  of  tidal  analysis  ^tplied  to 
discrefized  dau  given  at  (possibly  random)  times  t  >  1. 

K,  utilizes  (1)  in  the  form 

v(»*)  »  +«-  {1,2,  VP„*?.(t*-T.)}+  £(/*)  (3) 

where  the  summation  over  m  ranges  from  1  to  M  while  that 
over  n  ranges  from  I  to  N  and  where  the  following  defini¬ 
tions  are  employed.  (1)  The  set  (T^)  is  a  set  of  Af  reference 
series  constructed  from  the  tidal  constituents  obtained  either 
from  a  harmonic  analysis  of  a  lengthy  tidal  pressure  record  or 
from  tidal  potential  theory.  For  example,  coiutituent  ampli¬ 
tudes  Aj,  frequencies;^  and  phases  form  a  single  reference 
series  when  corniced  as  *  ljAjtxp{H2xfjt^-*fp], 

where  the  sum  in  J  ranges  ova  ite  constituenu  from  to  j, , 
which  forms  the  mth  of  the  linear  set  of  reference  series.  Su^ 
a  set  was  used  in  this  study.  But  equ^  accqNable,  w  need 
dictates,  are  reference  sets  fonned  by  individual  lines  or  by 
bilinear  combination  of  lines.  (2)  The  set  (t,)  is  a  set  of  N 
times  (including  zero)  by  which  the  reference  series  set  (T^) 
is  lagged.  (3)  Wq,  a  real  number  representing  the  mean  value 
of  the  series,  and  of  MN  complex  numbers,  are 

the  response  weights  that  are  a  measure  of  the  response  of  the 
velocity  series  v  to  the  reference  series  T^.  Equation  (3)  has 
2AfAf'4-l 'weights  (or.  equivalently,  degrees  of  freedom)  for  the 
fiL  (4)  The  series  e(fj^}  is  a  time  series  of  K  elemenu  con¬ 
taining  the  residual  error  between  the  dau  and  the  lit 

The  number  of  reference  series  and  the  number  of  lags  (Af 
and  Af,  respectively)  are  chosen  according  to  the  requirements 
of  the  problem.  For  the  analysis  presented  in  this  article,  all 
diurnal  corutituentt  are  summed  ^  all  semidiurnal  constitu¬ 
ents  are  summed  so  that  Af  •  Z  However.  Af  may  be  chosen 
suitably  large  either  for  a  treaunem  of  shallow-water  tidal 
effectt  analogous  to  extended  harmonic  analysis  [Zelfcr  and 
Cummings,  1967]  or  for  an  qrplicaiion  of  bilinear  analysis 
[Munk  and  Cartwright,  1966].  While  the  qrpropriate  num¬ 
ber  of  lime  lags  N  will  be  more  completely  addressed  in  the 
next  section,  note  that  the  fits  shown  in  this  article  utiUze 
three  lags  corresponding  to  -4S,  0,  and  448  hours  [Cart¬ 
wright  el  ai„  1969]. 

To  this  point,  the  selective  method  mimics  the  standard 
method  closely.  With  the  least  squares  minimization  of  the 
resktoalc.  however,  the  methods  ^erge.  For  the  selective 
method,  the  weights  Wq  and  IW^l  are  computed  from 

min  I^{S(i*)Ic(t*)]* ).  (4) 


where  the  suiiunation  over  k  ranges  from  1  to  AT  and  5  is  a 
real  series  defined  for  all  A  as 

5(1^)  >  1.  “accqMable'’ datum  at 

(5) 

5(1^)  *  0.  “tmaccepiable'‘datumaita 

which  weights  the  least  squares  fit  of  the  data  so  as  to  include 
{S  *  1)  or  exclude  (5  s  0)  a  datum  from  the  detomiiution  of 
the  response  weights.  The  weights  Wg  and  fit  only 

those  dau  selected  by  5  as  acceptable  and  is  the  mean  of 
only  those  data.  Minimization  (4)  is  entirely  equivalent  to  the 
least  squares  minimization  of  le(/j)l  over  the  A  >  1,  ...  K  ' 
dau  points  which  form  the  subset  of  accepuble  dau  defined  by 
identity  (S). 

Depending  upon  the  character  of  the  dau  series  to  be  aiu- 
lyzed.  construction  of  S  may  be  quiu  subjective,  demanding 
careful  attention  aitd  physical  intuition.  One  obvious  instance 
of  unacceptable  dau  occurs  with  sensor  malfunction  (e.g.. 
missing  dau  or  instrument  glitches).  A  less  obvious  but  very 
useful  irutance  might  be  the  contamination  of  a  velocity 
record  by  anomalous  evenu  (e.g.,  passage  of  a  bore)  whidi  the 
analyst  wishes  to  sqrarate  from  the  tidal  signal. 

The  minimization  required  by  (4)  can  be  performed  by  any 
standard  method  [Lawson  and  Hanson,  1974]  although 
some  are  numerically  superior  to  others.  In  particular,  the  nor¬ 
mal  form  of  the  least  squares  equations  may  result  in  numerical 
uttUbilities  that  can  cause  problenu  in  the  response  aiulytis 
[Adunk  and  Cartwright,  1966]. 

The  second  modification  of  the  standard  response  method  of 
tidal  analysis  involves  the  "opUihum  wigglinm*  criterion  of 
Zeller  and  Munk  (197S]  for  selection  of  N,  the  number  of 
lag  times.  Hie  criterion  rests  upon  the  variance  ratios  defi^ 
for  each  tidal  band  u 

r(A0  -  a*(E(A0)/a*(v).  (6) 

where  97(c)  am)  a^{y)  are,  respectively,  the  variance  in  the 
residual  e  of  the  fit  and  the  variance  in  the  dau  v.  Since  e  is 
a  fuciion  of  the  number  of  lags  or  weighu  N  employed  to  fit 
the  data,  r  also  dqiends  upon  N. 

Two  muniaDy  exclusive  pieces  of  dau  of  equal  length,  labeled 
A  and  B,  are  selected  from  the  record  to  be  analyzed.  For  a 
specified  number  of  lags  N,  the  weights  calculated  from  record 
A  are  enqrloyed  to  fit  record  A  and  to  predict  record  B.  The 
residuals  of  the  fit  and  of  the  prediction  yield  the  variance  ratio 
for  each,  namely,  and  respectively.  This 

procedure  is  repeated  using  the  weighU  from  dau  record  B  to 
obtain  the  ratio  r^^\N)  for  the  fit  of  B  and  r^^\N)  for 
the  predition  of  A.  These  variance  rations  may  be  computed 
separately  for  the  diurnal  and  semidiurnal  tidal  bands  by 
summing  the  squares  of  the  Fourier  coeffidenu  assodated  with 
frequencies  within  each  separate  band.  The  ratios  obtained  are 
ro^iNo).  roiAB^No).  ro^^KNo).  awI  rDVA)(S^) 
the  diurnal  band  and  rgiA^s),  rs^ABH/dg),  rgWQig), 
and  rgiBA)(fi/g)  for  the  semidiurnal  band. 

The  selective  method  must  be  slightly  modified  from  the 
procedure  just  outlined.  Because  some  dau  mi^  be  excluded  by 
the  dau  selection  series,  the  selective  response  method  cannot 
use  Fourier  analysis  to  calculate  variance  ratios  independently 
for  the  diurnal  semidiurnal  tidal  bands.  Rather,  the  vari¬ 
ance  ratios  must  be  computed  over  the  entire  frequency  spec¬ 
trum  from  the  fundamentd  definition 
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-  <«-  <£>f>  / <(v  -  •cv>)*>.  (7) 

who*  a  bndielMl  lime  (cries  is  simjdy  the  avenge  over  the 
IT*  aoceptaUe  data  points  of  the  series. 

For  both  die  standard  and  the  selective  medwds  of  tidal  anal* 
ysU.  it  is  expected  that  as  N  increases.  and  will 

decrease,  the  fiu  improvinc  with  an  increasing  number  of 
degrees  of  freedom.  Increasing  the  number  of  weights  utilized 
allows  mcee  aperiodic  noise  to  be  included  in  the  fit  On  the 
other  hand,  ratios  and  will  not  necessarily 

decrease  with  inaeasing  N  because  the  noise  of  one  record 
(which  determined  the  weighu)  is  not  predicuUe  in  the  other 
record.  According  to  the  optimum-wiggliness  criterion,  the 
number  of  diurnal  and  semidiurnal  weighu  and  Ng, 
respectively)  which  should  be  chosen  are  those  that  minimize 
the  variance  ratios  of  those  predictions  or  and 

or  respectively.  Similarly,  the  value  of  N  that 
should  be  chosen  for  the  srieedve  method  is  that  value  which 
minimizes  or 

The  essence  of  the  modificatioiu  to  the  standard-response 
method  of  tidal  analysis  may  be  summarized  as  follows.  (1) 
Instead  of  utilizing  all  dau  poinu  for  a  least  squares  Cl  only 
those  selected  by  S  as  acceptable  contribute  to  the  calculation 
of  the  response  weighu.  Consequently,  contamination  of  the 
fit  (including  the  mean)  by  evenu  and  missing  dau  can  be 
avoided.  (2)  Because  the  calculation  of  the  variance  ratio 
selecu  only  certain  data,  the  choice  of  the  number  of  lags  w 
respoiue  weighu  for  each  frequency  band  must  reflect  tidal 
behavior  across  the  entire  fri^uency  spectrum  rather  than 
being  determined  separately  for  each  of  the  tidal  bands. 

EXAMPLES 

Attention  is  now  turned  to  a  comparison  of  the  standard  and 
selective  response  methods  of  tidal  analysis  when  each  is 
applied  to  the  daU  set  that  initially  prompted  the  modifi¬ 
cations  outlined  in  section  2.  These  Doppler  curreni  records 
comprise  three-dimensional  velociqr  measurements  with  a  10- 
m  vertical  averaging  interval  over  a  depth  range  of  40  m  to 
240  m  and  with  a  S-min  time  averaging  iniervaL  The  observa¬ 
tions  were  made  in  die  Tarifa  Narrows  of  the  Strait  of  Gibraltar 
from  March  17  through  April  23,  1986. 

All  of  the  analyses  presented  in  this  section  use  a  set  of  ref¬ 
erence  series  derived  from  a  long,  coastal  tidal-elevation  su- 
tion  (Tarifa,  Spain)  record  and  a  6-tnonih  time  series  of  bottom 
pressure,  fri  e^  instance,  three  time  lags  of  two  summed  tidal 
bands  (N  >  3  and  Mm2)  provide  13  (2AfAr.fl)  degrees  of 
freedom  for  the  least  squares  weighu  (including  the  mean). 

The  first  exam{de  is  the  attempt  of  the  standard  and  the  selec¬ 
tive  response  methods  to  extract  the  tide  and  mean  from  a  very 
problematic  time  series  of  currenu  conuining  both  dau  gqis 
and  large  eventt  associated  with  an  inumal  bore.  Figure  1 
dis{days  a  brief  portion  of  the  Doppler  dau  fin  the  eastward 
component  of  velocity  ((/)  at  a  dq>th  of  40  m.  This  record  is 
shown  in  Figure  Id. 

In  this  particular  dau  record,  it  is  difficult  visually  to  distin¬ 
guish  the  intenia)  bore  from  the  tidal  oscillation.  However, 
the  signature  of  the  bore  in  other  daU  is  well  defined  (Pef- 
ifgrcw  and  Hyde,  1990].  For  example,  at  dq^  new  the 
pyc-nocline,  iu  velocity  front  is  sufficiently  sh^  itt  arrival 
is  sufficient^  jdiase  delayed,  and  itt  duration  is  suffi^ienlly 
brief  that  it  is  easily  differentiated  from  the  tide.  Moreover, 
die  displacement  of  the  pycnocline  fay  the  passage  of  the  bore 
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Fig.  1.  CompadMe  of  the  nandaid  aid  selective  letponse  meihods  of 
tidal  analyat  as  spplied  lo  Tarifa  Ooppier  dau  (esneraid  velocity.  40- 
as  depth)  widi  eveau  sod  gaps.  Velodiy  onitt  aic  cm  s'>.  (a)  Velocity 
time  series  with  the  gaps  filled  by  imeipolation  (solid  line)  and  the 
aiean.phis.tide  fit  by  the  standard  method  (dotted);  (b)  lesidiial  by 
dw  standard  method;  (c)  dau  selection  series  for  the  selective 
method;  (d)  velocity  time  series  (solid  line)  and  the  mean.p)ns4ide  fit 
by  the  sde^e  method  (doited  line);  (e)  icsidual  by  the  selective 
methodL 


aids  m  delineating  die  surge  of  the  bore  from  the  tidal  cutrenu. 
Such  corroborating  evidence  indicates  that  the  sudden  rise  in 
velociqf  from  roughly  80  cm  s'*  lo  over  ISO  cm  s'*  during 
each  tidal  cycle  is  associated  with  the  velocity  from  of  the 
inlemal  bore.  At  40-m  depth  in  the  Tarifa  Narrows,  the  arrival 
of  the  front  occurs  slightly  after  the  maximum  (east-ward)  flow 
of  the  underlying  tidal  oscillation.  During  the  pas-sage  of  the 
bore,  the  velocity  remains  elevated  for  almost  6  hours  before 
falling  rapidly  to  values  of  roughly  70  cm  s'*.  Following  the 
bore,  there  is  usually  a  period  of  missing  data.  These  gqw  are 
associated  with  vertical  excursions  of  die  instrument  (mounted 
on  the  upper  flotation  sphere  of  a  subsur-face  mooring)  to 
dqMhs  from  which  sampling  at  40  m  was  unattainable. 

The  daU  are  first  fitted  using  the  standard  response  method. 
The  relatively  short  daU  gaps  we  filled  unobtrusively  using  a 
linew  interpolation.  A  mean  and  weighu  we  calculated  using 
the  entire  dau  record  including  the  dau  wsociated  with  the 
internal  bore.  Figure  lo  diqilays  the  fit  of  the  mean  plus  tide, 
and  Figure  16  shows  the  residual  of  the  fit  While  the  fit 
qipears  lo  be  good,  it  includes  not  just  the  tidal  oscillation  but 
also  the  motion  identified  with  the  surge  of  the  internal  bore. 

The  selective  response  method  of  tidal  analysis  provides  a 
very  different  fit  If  the  analyst  considers  the  dau  associated 
with  both  the  velocity  surges  (pwsages  of  the  bore)  and  the 
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dau  gaps  (divas  of  tlia  inamimont)  lo  be  unaccapiabic,  ihaa  lha 
satectioa  tarias  5  it  constructed  at  in  Figuia  Ic.  The  fit  geii- 
eralad  by  dia  telacdva  rasponsa  method  (which  uses  only  dau 
pointt  conasponding  to  the  conact  level  of  the  instnimcai  and 
to  the  absaaca  of  the  bore)  it  displayed  in  Rgure  Id.  Figure 
la  it  the  residual  of  d  li  analysis. 

Application  of  the  two  methods  has.  not  too  surprisingly, 
led  to  drastically  dilfmat  interpreutions.  Compare  Hgure  lo 
to  ingura  14  and  Figure  Ih  to  Figure  la.  The  result  of  the 
standard  requose  mediod  suggetu  a  tidal  anqjditude  of  roughly 
7S  cm  t**  superimposed  upon  a  mean  flow  of  ioughly  SS  cm  s'* 
with  a  reridual  dominated  by  tidal  harmonict  having  a  mean  of 
zero  nd  an  amiditude  of  roughly  23  cm  s'*.  In  contrast,  tha 
selective  response  method  yields  a  smaller  tidal  oscillation 
(roughly  30  cm  s*’).  a  weaker  steady  flow  (roughly  30  cm  t*'). 
aitd  an  enormous  potitive.definiu  residual  (roughly  33  cm  s'* 
average  cantriburimi  to  the  mean  flow).  Moreover,  tha  fits  by 
the  two  methods  show  a  phase  diflcrence  of  nearly  3(F  with  the 
fit  of  the  standard  responm  method  shifted  toward  the  advent  of 
the  bora. 

Supporting  evidence  (for  example,  analyses  at  depths  where 
the  bon  is  unmistakably  defirud)  clearly  indicates  that  in  this 
particular  case  the  picture  presented  by  the  selective>response 
method  is  more  indicative  ot  the  local  dynamics  than  that  pre¬ 
sented  by  the  standard  method  [see  Petiigrtw  and  Hyde. 
1990}.  The  standard  method  fails  lo  portray  the  positiva  defi- 
nitt  luoira  of  tha  events,  erroneous  tidal  harmonics  appear  as 
an  artifact  of  the  aiulysis,  and  the  surges  of  the  mlemal  bora 
contaminau  the  value  of  the  underlying  steady  flow.  By  con¬ 
trast.  the  selective-response  method  avoids  these  pitfalla. 

Tha  second  example  tasu  tha  alnltty  of  tha  aelocriva 
response  method  of  ti^  analysts  lo  meet  the  optifflum-wig- 
gliness  criterion  for  selecting  the  sppropriau  number  of 
weights.  For  the  test,  two  mutually  exclusive  series  of  34S- 
hour  duration  each  (approximately  a  single  spring-neap  tidal 
cycle)  are  selected  from  dau  of  the  eastward  component  of 
velodiy  at  a  dqidi  of  240  m  in  Tarifa  Narrows.  The  series  are 
denoted  A  and  B.  respectively.  At  240  m.  the  influence  of  the 
bora  is  relatively  minor  and  there  are  no  gape  in  the  data.  The 
flu  and  residi^  by  the  two  methods  consequently  have 
relatively  slight  differences.  Thus  this  exam^  provides  a 
reasonable  comparison  of  the  optimum  number  of  weights 
selected  fay  each  mediod. 

As  in  the  previous  example,  die  selective  reqioiise  method 
calculates  a  inean-plus-tide  fit,  rejecting  daU  when  the  bore  is 
known  to  be  present;  the  dau  selection  series  iherelbre  resem- 
tdes  the  curve  in  Rgure  le.  The  standard  method  uses  aO  data. 

Recall  that  the  general  criterion  for  optimizatioa  is  based 
upon  the  ratio  of  the  variance  in  the  residual  to  the  variance  in 
ite  dau  signal  m  defined  in  (6).  The  performance  of  die  selec¬ 
tive  response  mediod  is  periiqie  most  fairly  assened  by  com¬ 
paring  the  variance  ratio  calculated  by  the  selective  method 
''gfN)  (which  encompasses  the  entire  frequency  qiectmm) 
with  the  variance  ratio  obtained  from  the  standard  method 
'ijgQf)  when  the  Fourier  componems  of  both  the  diurnal  and 
the  semidiuinal  bands  are  summed  togetlm. 

For  the  present  examjde,  the  length  of  the  dau  records  results 
in  the  use  of  26  (conqrlex)  Fourier  components  for  the  com¬ 
bined  diurnal  Odl  ±  4  J  cpm  or  05664  ±  0.1647  epd)  and 
semidiumal  (JMl  1  45  cpm  or  1.9323  ±  0.1647  cpd)  bands. 
The  energy  of  the  semidiumal  tide  oonqiletely  donrinates  diat 
of  die  dhrmal  tide  so  that  m 

Figure  2  dwws  that  the  two  methods  have  produced  eqniva- 
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Fig.  2.  Caesparison  «f  Ike  optiniuw  wiggliaeM  ciiieriee  u  datci- 
auaed  by  Ike  tuodaid  md  relive  icspeotc  wrtkoJi  applied  le  two 
ccnaecodve  tpaag-eeap  cycles  (A  and  B).  RadM  aie  ef  Ike  variaaoe 
ia  Ike  leaidnal  to  die  variaoee  is  ike  ai|^  (6)  sa  a  liuctkiB  ef  ike 
aatebat  ef  weighis-  Rsimt  fat  ike  awnaeed  dkmal  and  senidiaBaal  6c- 
qsesey  baoda,  which  were  obieined  fiem  cycle  A  weigbu,  arc 
(fir  ef  A)  and  rMteto  (pndicricai  ef  Bk  iboae  bam  cyde  B  weaghu  arc 
rmtas)  (find  a)  mi  rp/tM  (pmSmboa  ef  A).  Ratios  for  the  catire 
ftaqetnry  apaeuwn.  which  wen  obiaintid  from  cycle  A  weights,  are 
rg(M)  Qh  of  A)  and  rnUS)  (pradictioa  of  Ak  ibose  from  cycle  B 
woi|ktt  are  r  ^lam  ^of  B)  wd  (pmfiotice  of  A). 


lent  results,  the  dependence  of  on  N  being  roughly  that  of 
rgf  on  N.  Specifically,  the  self-fit'  variance  ratios  (namely, 
ridA  and  deen^  with  increasing  N.  This  improve¬ 
ment  of  the  fit  is  of  course  to  be  expected  with  an  increasing 
nu«">wr  of  degrees  of  freedom.  However,  the  predicted-fit  vari¬ 
ance  ratios  (namely,  r4AJ)  r(*^>)  show  an  eventual  dete¬ 
rioration  of  the  fit  with  an  increasing  number  of  response 
weightt.  So,  following  the  ptescripiioa  of  Zeller  and  Munk 
[1973]  for  this  examide,  both  methods  suggest  that  the  opti¬ 
mum  number  of  weightt  is  two  for  the  AB  ratios  (that  is,  the 
prediction  of  cycle  B  from  the  weightt  of  cycle  A)  and  frar  for 
the  BA  ratios.  A  compromise  of  three  lag  times  (Af  ■  3)  would 
seem  reasonaUe  for  this  case. 

It  can  be  concluded  that  use  of  the  ratios  rg(^k)(N)  and 
of  the  selective  response  method  would  produce 
resultt  reasonably  equivalent  to  those  of  the  standvd  response 
method  when  the  variance  outside  the  diurnal  and  tenudiumal 
frequency  bands  is  small,  te..  when  r^CN)  and  roy(N)  are 
qniu  similar.  Becaufe  the  selective  method  must  utilize 
r|j(A0  for  this  analysis,  it  is  unable  to  determine  a  diflerem 
optimal  number  of  weightt  for  the  various  sqiaraie  tidal  bands. 

DBCUSSXm 


The  selective-response  method  of  tidal  analysu  employs  a 
slight  modification  of  the  standard  method  to  inqvove  consid- 
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erably  the  extraction  of  the  steady  and  tidal  components  from 
proMematk  velocity  time  series.  A  dau  selection  series  (5  in 
definition  (S))  simply  zero  wei|hts  missing  or  anomalous  dau 
in  the  least  squares  determination  of  the  response  weights. 
Consequently,  the  tidal  analysis  remains  uncontamiiuted  by 
such  features. 

The  selective  response  method  shows  promise  in  dealing 
with  dau  gaps.  Because  the  fit  U  not  forced  to  conform  to  an 
artificial  interpolation  of  the  missing  data,  the  globally  dam¬ 
aging  effects  of  interpolation  are  never  introduced.  Rather,  by 
ignoring  gaps,  the  method  permits  the  response  weights  and 
the  reference  series  to  extertd  the  fit  naturally  across  them. 
Numerous  test  cases  (not  presented  in  this  article)  indicate  that 
the  selective  method  is  capable  of  handling  missing  daU  not 
only  for  long  periods  (many  tidal  cycles)  but  also  for  frequent 
loues  (more  than  half  of  the  record).  Quite  possibly,  the 
method  proposed  in  this  article  will  find  its  widest  aitd  most 
useful  application  to  be  that  of  bridging  dau  gaps. 

The  selective  response  method  can  also  be  useful  in  isolating 
anomalies  or  evenu  in  a  current  record  that  may  be  phase 
locked  to  the  tides.  As  the  example  of  the  internal  bore  in  the 
Strait  of  Gibraltar  demonstrated,  the  method  also  improves  the 
determination  of  the  steady  and  tidal  flow  components.  Spe¬ 
cifically,  the  inclusion  of  dau  associated  with  non-zero- 
averaged,  phase-locked  surges  (such  as  those  of  an  internal 
bore)  may  force  the  standard  method  to  estimau  poorly  the 
steady  flow  and  tidal  amplitude,  shift  the  phase  of  the  tide,  and 
leave  unresolved  the  structure  of  surges.  This  should  be  of  seri¬ 
ous  concern  in  standard  tidal  analyses  of  the  Strait  of  Gibraltar 
and  similar  environments  where  the  internal  bore  may  domi- 
nau  velocity  time  series. 

nnally,  it  must  be  noted  that  the  selective  response  method 
has  its  limitations  and  its  successful  usage  cannot  be  guaran¬ 
teed.  The  following  caveats  are  offered.  (1)  The  placemern  of 
ihe  unaccqitable  dau  is  a  critical  factor.  For  example,  if  the 
dau  gaps  oi  anomalous  events  occur  such  that  the  dau  selec¬ 
tion  series  zero  weights  dau  at  the  maxima  and/or  the  minima 
of  most  of  the  tidal  cycles,  then  the  method  will  not  be  able  to 
discern  the  correct  amplitude  for  the  fit  Test  cases  of  this  sort 
(not  presei.ted'here)  rarely  proved  reliable.  (2)  Objectives  of 
the  aiulysis  must  be  clear,  and  definitions  of  "mean"  and  “tide* 
must  be  precise.  For  instance,  in  the  analysis  presented  for  the 
Strait  of  Gibraltar,  the  expressed  goal  is  the  recovery  (as  a 
residual)  of  the  internal  bore.  The  "mean*  and  the  "tide*  there¬ 
fore  have  to  be  appropriately  defined  as  the  background  upon 
which  the  surges  of  the  internal  bore  are  superimposed.  (3) 
Anomalous  baroclinic  events  may  be  clearly  discemible  at 
some  levels  of  the  water  column  and  masked  at  othen.  Isolated 
observations  within  the  water  column  may  or  may  not  reveal 
these  features.  Thus  the  ability  to  construct  an  effective  selec¬ 
tion  tinie  series  may  depend  upon  good  fortuiw. 
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